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1. Introduction 

The strong interaction energy shifts ei s and to- 
tal decay width r ls in pionic hydrogen have been 
measured to a remarkable precision pp 

e u = [-7.108±0.013(stat)±0.034(syst)] cV, (1) 
r ls = [0.868 ± 0.040(stat) ± 0.038(syst)] eV. (2) 

It is well known |2I3| that the (complex) strong 
interaction shift in the Is state of hadronic atoms 
is closely linked to the (complex) elastic threshold 
scattering amplitude a defined in the absence of 
the Coulomb field. We refer in the following to 
this quantity as the (complex) scattering length. 
This is conventionally expressed in the ratio of 
the shift to the Bohr energy Eb = —ma 2 /2: 

|l = ^(l + y=4mQa(l + i l8 ), (3) 

where <5i s conveniently measures the deviation of 
the shift from the lowest order estimate 

4 --£««».. w 

Here <^s(r) is the non-relativistic Is Bohr wave 
function of a point charge and m the reduced 
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mass, which in the present case is that of the 7r~p 
system. It is important to understand the correc- 
tion Sis transparently and reliably to an accuracy 
matching the high experimental precision, since 
the hadronic 7rN scattering lengths are key test- 
ing quantities for chiral physics. In addition, they 
are needed phenomenologically to about 1 % for 
the precision determination of the 7rNN coupling 
constant using the GMO relation 4 . 

The standard conversion of experimental data 
to a scattering length uses the potential approach 
of Sigg et al. 5 , which describes the 7rN interac- 
tion in terms of coupled equations using phys- 
ical pion masses and an isospin invariant non- 
diagonal potential matched to scattering lengths 
calculated by setting the neutral pion mass equal 
to the charged one. This gives #i s (Sigg) = 
(— 2.1±0.5)%. The procedure is model dependent 
and it is not consistent with the 7rN low-energy 
expansion 4. Their results must therefore be 
used with caution. 

The classical way to obtain Eq. © is based on 
analytical approaches using Coulomb wave func- 
tions (see Refs. |2l3lb17l8] and references therein). 
To our knowledge these papers do not explore the 
effect of the extended charge distribution on the 
strong interaction shift. This paper discusses this 
question. 
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The authors of Refs. |9ll(Jj have calculated the 
ground state energy of the ir~p system in the 
framework of QCD+QED, using effective field 
theory (EFT) techniques. The shift of the ground 
state energy is related to the scattering lengths 
in pure QCD, evaluated in the isospin symmet- 
ric limit m u = md The corresponding correction 
5 e in the energy shift is evaluated in the frame- 
work of Chiral Perturbation Theory (ChPT). At 
leading order 0, one has S e =(-4.3 ± 2.8)%, 
while the next-to-leading order result JUj is 5 e — 
(— 7.2±2.9)%. The uncertainty in S e is due to the 
poor knowledge of one of the low-energy constants 
occurring in the effective theory. The correction 
8 € is also considerably larger than what was found 
in Ref. This EFT approach makes inherently 
no distinction between the atomic corrections due 
to Coulomb effects discussed here and other con- 
tributions. We come back in section 4 to a com- 
parison with the present work. 

Our aim here concerns only the connection of 
the strong atomic energy shift to the scattering 
length a h defined as the one which would be ob- 
served if the Coulomb field of the extended charge 
could be removed and considered as due to an ex- 
ternal source. This scattering length is directly 
related to the one appearing in forward disper- 
sion relations for 7rN scattering. In this spirit, no 
correction is made for the internal e. m. contri- 
butions to the masses. The physical scattering 
lengths for the 7r ± n scattering correspond closely 
to the present definition, neglecting the very small 
e. m. correction from the charged pion interac- 
tion with the neutron charge distribution. These 
scattering lengths with physical masses are natu- 
ral 'observables' for the study of isospin breaking. 
The (complex) 7r + n scattering length coincides 
with the 7r~p one in the limit of exact charge sym- 
metry: it has the corresponding open charge ex- 
change channel 7r + n — > 7r°p and the correspond- 
ing open radiative decay channel 7r + n — * 7p. This 
definition is different from that of the QCD scat- 
tering length used in the EFT approach. We use 
here the r-rcprcscntation, which is more transpar- 
ent for the present problem than the equivalent 
momentum representation. Since the 7r~p atom 
is highly non-relativistic, most of the discussion 
will be made using non-relativistic concepts. The 



result will be expressed in terms of the empirical 
on-shell parameters of the 7rN low-energy expan- 
sion. For the electromagnetic corrections, this 
approach gives intuitively interpretable expres- 
sions, exact up to terms in a 2 log a, provided the 
Coulomb potential of the extended charge varies 
little over the range of the strong interaction. 

In section 2 we solve the problem exactly to all 
orders in a in the limit of a short ranged strong 
interaction with the charge located to a spheri- 
cal shell. The correction for an arbitrary charge 
distribution is then derived perturbatively to the 
same order in a as in the EFT expansion |l()j . 
Corrections for the finite interaction range are ex- 
plored. We explicitly include the correction for 
the vacuum polarization. In section 3 we discuss 
the magnitudes of the corrections and their physi- 
cal structure. In section 4, we compare our results 
to those of previous approaches. 

We denote by Ej s ,E and e\ s — E — Ef s the 
Is finite size e.m. binding energy, the total Is 
binding energy with strong interaction and finite 
size, and the strong interaction shift, respectively. 
The non-relativistic wave numbers are Kg, Kf s 
and k for EB,Ef s and E, respectively. Since this 
paper concerns atomic corrections we also use the 
Bohr radius tb = >^b = ( ma ) 1 - 

2. A Model for the 7r p Atom 

The aim of this section is to explore the con- 
sequences of the extended charge distribution in 
a pedagogically transparent and soluble model of 
the 7r - p system. This will serve as a prototype for 
the later more general discussion and it will reveal 
the nature of the contributions to the correction 
term in Eq. J3J. In the absence of the Coulomb 
potential the threshold expansion for the angu- 
lar momentum I = 0, typical of a weak scattering 
length, is related to the phase shift Si = o and to 
the momentum q by the relation 



Here a h is the hadronic scattering length and b h 
is the range parameter. 

The model is constructed as follows. We first 
consider the case of a single channel. This avoids 
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the complications of several open channels with 
different masses for which the equivalent of the 
single channel scattering length is ambiguous. 
The generalization will be discussed later. The 
charge is assumed to be concentrated to a spher- 
ical shell of radius R, outside the range of the 
hadronic interaction. The system is taken to be 
non-relativistic. For the moment we neglect the 
effect of the vacuum polarization potential. As 
defined, this problem can be solved exactly in 
terms of the on-shell hadronic s-wave scattering 
amplitude, although we will only evaluate it for 
contributions to the correction term up to order 
a 2 log a. In this case the definition of the scat- 
tering length is clear. The model does not have 
problems with the intermingling of the Coulomb 
and hadronic interaction, contrary to other de- 
scriptions. Such problems are particularly acute 
in any description with a pointlike charge distri- 
bution, since the Coulomb interaction is then di- 
vergent at r = 0. 

Inside a typical shell radius R of the order of 
1 fm, the Coulomb potential is constant with 
V C (R) = -a/R ~ 1.4 MeV. This motivates the 
non-relativistic approximation in this region to 
order a. The inside wave number q c is constant: 



<17 



2ma 
R 



2ma 
R 



(6) 



The Is binding energy E ~ — 3.2 keV is negligible 
compared to the Coulomb field and the strong 
interaction inside the charge distribution region, 
although its exact value governs the scale of the 
atom. The external Is wave function for r > R 
is a Whittaker function Wa ; i/2( z ) ( see e - g-> |TT) . 
Eq. (9.237)) with A = k b /k and z = 2nr. The 
inside wave function for r < R is a standing wave 
outside the strong interaction region. Neglecting 



terms of order(l — A) 2 ~ (ama h ) 2 
has for r < R, 



u ln (r) = N 



sin (q c r) tan.(S£_ n . . 
+ cos (q c r) 



(7) 



and for r > R, 

u out {r) = (4 7 r) 1 / 2 )- 1 Kexp(-z/2) 



x{4(l+(l-A)(l- 7 -log 



(8) 



Here 7 = 0.577... is the Euler constant. Note that 
the term tan(5 i h =0 ) is determined by the 'phys- 
ical' hadronic phase shift in the absence of the 
Coulomb field taken at the energy —Vc{R)- The 
wave function corresponding to Eqs. (7J and JHJl 
is normalized to order a 2 . 

The energy shift t\ a produced by the strong in- 
teraction is obtained by matching the logarithmic 
derivative of the wave function at the radius R. In 
accordance with standard practice, it is defined as 
the difference between the total binding energy E 
and the electromagnetic binding energy with a fi- 
nite size charge distribution |S] . This corresponds 
to the removal of the 'scattering length' corre- 
sponding to the extended charge distribution: 



amR 2 
a fs = — ; e 



(9) 



Expanding the exact analytical expression to 
terms of order a 2 log a by a straightforward al- 
gebraic calculation gives the following correction 
factors, where the hadronic scattering length a h 
takes the place of a in Eq. J3J: 



Sis 



n R n a h r„ 
-2 \-2 — 2- 

Vb r B 
2ma b h 



7 — \og(2amR)] 



R a h 



(10) 



This expression serves as a guide for the later gen- 
eralizations. 

The assumption of a zero range hadronic inter- 
action is unnecessary in our simple model. An 
interaction of any shape will give the same result 
provided its range is smaller than R. This follows 
from the matching condition for the wave func- 
tions and ||BJ|, which is only required at R, 
such that any interaction with the same hadronic 
scattering near-threshold amplitude, a h + q 2 b h , 
gives the same result. 

The terms in Eq. i|10[) have a clear physical in- 
terpretation. The extended charge wave function 
at r = in the absence of strong interactions 
is (j>in(0) = 0s(O)(l — R/tb + ■■) to the present 
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order in a. It is a better starting approxima- 
tion than the wave function of the Bohr atom 
in Eq. J3J, which then should be multiplied by a 
factor (1 - R/r B + -f ^ l-2R/r B + - This ac- 
counts for the first term in Eq. (|10|l . The second 
term proportional to a h is a renormalization due 
to the external wave function which is changed 
at R by the hadronic scattering itself by a factor 
l + 2a h ma[2 — 7 — log(2ami?)]. The outside wave 
function is determined by the energy shift. The 
matching of the inside and outside wave functions 
of Eqs. J7J) and (JSJ gives near the origin 

U in (r) = {l + 2a /l ma[2-7-log(2ami?)]} 



x (r + a h ) m (O). 



(11) 



This result agrees with that obtained by match- 
ing the logarithmic derivative at R. This factor 
has little sensitivity to the exact value of the ra- 
dius R. The leading a log a part of the term in 
a h is well known from previous approaches and 
has also been found in the EFT approach where 
it corresponds to a "loop" term |9ll()j . 

The last term in Eq. (|1C)|> follows from gauge 
invariance with the replacement E — > E — Vc(0) 
in the scattering amplitude g]. Alternatively, 
and more intuitively, it follows using the correct 
energy at the point of interaction. This is not 
the binding energy, but the finite depth of the 
Coulomb potential of the extended charge. (For 
the corresponding effect in higher Z pionic atoms, 

see Refs. [12- D3). 

Exactly the same reasoning as for the Is state 
can be applied to hadronic energy shift e ns in 
states of any n. The correction factor 5„ s is de- 
fined in complete analogy to Eq. Q: 



Ens 



E,, 



"(I 



4 (l + 6„ 

nr B 



(12) 



where E ns = —ma 2 /(2n 2 ) and the convenient 
comparison shift e ns is the simplest perturbative 
expression for the energy shift 



2m 

One has to order a 2 log a in the correction 



= 2— M--+5J log n 
r B \ n ^ k 



(13) 



(14) 



In this expression all the dependence on the pa- 
rameter R has disappeared, which reflects that 
all the short-ranged physics is identical but for a 
renormalization factor. In the limit n — > 00 



5i s 



-(1 + 7): 



(15) 



where the correction term is given by 



-2 h 2 — 

TB TB 

2ma b h 



1 - 27 - log 



2R 



I'B 



R a 



h ' 



(16) 



This semi-classical limit for /t„ e (kb/b) ^ 
corresponds to the Coulomb scattering length ^B] 
a c = a h (l + <5 00s) in our model. For the present 
case of a ir~p atom the numerical difference in 
the correction terms for different values of n in 
Eqs. {TJ} and {T51) is less than 10 -3 and of little 
practical importance. 
Arbitrary Charge Distribution 
The result i|l(J|l is the prototype for more gen- 
eral charge distributions. The difference between 
the Coulomb potential Vcf s for a charge distri- 
bution from the observed n~ and proton form 
factors, p(r), and Vcr corresponding to that for 
the spherical shell of radius R, pn{r), gives a per- 
turbative potential, which includes the f.s. charge 
density associated with the anomalous magnetic 
moment: 

SV c (r)=Vcf s (r)-V CR (r) 



OO 

<j (-J--^j Sp(r'Hirr' 2 dr', 



(17) 



where Sp(r') = p(r) — pR(r'). Applying this per- 
turbation to our soluble model gives a net cor- 
rection independent of R to the present order in 
a. The explicit calculation leads to the follow- 
ing four changes in our model results. First, the 
e. m. finite size energy shift © is changed with 
the substitution of the model R 2 by (r 2 ) = 

J \ / em 

( r l)em + ( r ")e m = L15 ( 2 ) ^ &S in Ref ' ^ 

Likewise, the value of the overall Coulomb poten- 
tial at the origin changes from the model value 
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a/R to a (l/r) em and the wave function squared 
at the origin changes it value from (1 — 2aR + ..) 
to (1 — 2a (r) + ••)• Finally, the term logi? is 
replaced by (log r) em . The changes are indepen- 
dent of the hadronic interaction. 

In the case of a single channel and in the 
hadronic zero range limit, the corrections are: 



Sis = -2- 



2 ^( 2 _ 7 _/ log * 

TB \ \ T B 



+2ma 



A\ h = S {r) +5^ + 33. (18) 

Warn a 



We now introduce the correction 6 vp for the vac- 
uum polarization 16 . The first order vacuum 
polarization is described by a potential propor- 
tional to a 2 with a range much larger than that 
of the hadronic and charge distribution ones and 
it is insensitive to the strong interaction dynam- 
ics. The joint extended Coulomb potential and 
vacuum polarization one is a perfectly justifiable 
alternative to the point Coulomb potential as the 
starting point for the wave function in Eq. ©. 
The square of the unperturbed wave function at 
the origin changes by S vp — 28(f) vac {Q)/ 0b (0) = 
0.48% due to vacuum polarization J^j and by 
—0.85 % from the extended charge (see Tabic 0). 
This result is model-independent and it agrees 
with the prior numerical value implicit in Ref. [JJ]. 

In the derivation of the correction factors, we 
nowhere used that these quantities should be real. 
We can therefore take the energy shift to be com- 
plex, 6\ s — iri s /2, with a hadronic complex scat- 
tering length + ia^. The complex energy shift 
is related to the corresponding correction factors 
by a(l + S) — > a r (l + 6 r ) + itti(l + Si), as in 
Eq. (J3J). The imaginary part corresponds to ab- 
sorptive phenomena. In the notation of Eq. I|18|) : 



ois,i — —2 



<r> 



tb 



2ma ( — 
\r 



r ,h 



7 



log- 



2r 

TB 



5 {r) +2S c + S B i . (19) 



Here the imaginary amplitudes and b\ refer to 
any absorptive channel such as the 7r~p charge 



exchange scattering. Note the additional factor 2 
in the middle term as compared to that for the 
real case in Eq. (|18f) . Since a.j << \a r \, the change 
in 5i S)r due to absorption is negligible. 

We conclude that most of the corrections to the 
width are due to the change of the wave function 
at origin: it is important to use wave functions 
corresponding to the finite size and vacuum po- 
larization potentials. In addition, the non-linear 
renormalization term must also be included, but 
only the real part of the scattering length is rel- 
evant. To these should be added the amplitude 
change due to the gauge term in analogy to the 
case for the energy shift. 
Coupled channels 

The Tr~p atom is a coupled system of the con- 
tinuum ir°n and "fn channels in addition to the 
7r~p one. These three channels are denoted by 
indices i (J) — c,o,f, respectively. The low- 
energy expansion in multiple channel systems is 
defined in terms of energy dependent (symmetric) 
K-matrices which enter the standing wave solu- 
tions. The formalism is described briefly below 
and it is illustrated for the 2-channel situation. 
The single channel becomes a special case. The 
standing waves at distances larger than the charge 
radius r > R are defined as ^7] 



(20) 



where Kf 4 are the "Coulomb-corrected" K- 
matrix elements. The wave functions u c r and u c s 
are defined in terms of the standard regular and 
singular Coulomb functions F and G, respectively 
(see e.g. Refs. [till 7) ). In the limit a — * and q 
fixed, these solutions correspond to sin(qr)/g and 
cos(qr)/q, respectively. Furthermore 



u c _ F ( r ) _^ sin M . 
r C(rj)q r^oo C(r/)q' 

u c s = G(r)C(rj) - 2 V h( V ) 
-> C{rf) cos(<p) - 



(21) 



F{r) 

c(v) 

2r/h(r]) sin(<^) 



CM ■ (22) 

where r/ = zz'am/q and qrj = —kb- In 
these equations the digamma function ip(x) de- 
fines h(rj) = [^(irj) + ip(—ir))]/2 — log(?7 2 )/2 and 
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C 2 {rj) — 27r?7/[exp(27r?7) — 1] is the standard pen- 
etration factor. At large distances the phase 
if = qr — r]log(2qr) + a, where er is the Coulomb 
phase shift. 

From the iT-matrix one obtains the scattering 
amplitude T by regrouping the standing waves 
into the regular and outgoing waves. With C c = 
C(rj) and C Q — Cf = 1, one has 



1,3 



Ci[-{K c r i + f] 1 Cj 



(23) 



In this equation we use a diagonal matrix with 
fee = 2rjq c h{r]) + iq c C 2 (rj) and / ,o = iqo- 

In the single channel case, the textbook relation 
of K c to the scattering amplitude T is JS] 



r c = c 2 (r,) 



i 

~K~ C 



2rjqh(ri) + iqC 2 (rj) 



= — [qcot(S) — iq]~ 



(24) 



The atomic level shift is obtained from the "well 
known" formula. 



C 2 {rf)qcot{5) + 2rjqh(r]) 



1 

K~ c 



and the bound state condition cot(5) = i 
found by Trueman [3] 



47T 

2m 



0l(O)a c 



1 



-(2 + 27) 



(25) 
As 

(26) 



where a c is given by the Coulomb K-matrix at 
threshold. Its relation to the hadronic scattering 
length of the present model is 



1-- 



2R 



-a' 1 A c 



Ac 



2 



1 - 27 - log 



2ma 
2R 

TB 



(27) 
(28) 



XXX This relation leads to the correction given 
by Eq. JTSJ. 

In the 2-channel case with = (c, o), the 

leading order in the level shift follows by the re- 
placement of a c by the threshold amplitude in 
the charged channel A c c c , which is obtained from 
Eq. (|23|1 . Including terms to order q 2 



At the charged threshold the phase space left in 
the open neutral channel is described by the mo- 
mentum q a . Eq. (|26|l should now be used with a 
complex A c c c : 



Ait 
2m 



-iT ls /2 



1 



A': 



>B 



27) 



(30) 



In the zero range limit, the hadronic interaction 
in the charged channel occurs at a momentum 
q 2 = 2ma/R, while that in the neutral channel 
still occurs at the momentum q Q , since the atomic 
binding energy is negligible. The K-matrix el- 
ements are energy dependent with a low-energy 
expansion is analogous to that of Eq. (J5J . 



K c , c 

K ,o 



i n +6" 



2ma 
IT' 



< + ^{qI + qD bio, 



(31) 

(32) 
(33) 



We assume isospin invariance for the range pa- 
rameters 6j, l c and b^ since they only appear in 
correction terms. In the internal region r < R 
the standing waves are 



u) tx sin(qir)6ij + K hj cos(qjr). 



(34) 



The continuity of the wave function matrix it and 
its logarithmic derivative u~ 1 du/dr at the radius 
of the charged shell R gives: 



K c 



1 



2i? 

— + Kc, c A G 
rs 

R 

TB 



= K n 



K , c A G K C 



(35) 

(36) 
(37) 



These corrections are implicit in the single chan- 
nel equations (fTUfl and i|TT)|) . 

The extension to the (7, n) channel is ob- 
tained with the substitution if? ■ — > Kf j + 
(iqfK?jKf j)/(l — iqfKf j) to every matrix el- 
ement of the two channel system. The higher 
order terms in the neutral and photon channels 
K C Q and j are negligible such that 



A% c = Kl c + iq {Kl ) 2 - q 2 K c jKl ) 2 



(29) 



K c 

C,C 



+ zq {Kl ) 2 +iq f {Kl f Y 



(38) 
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The corrections (|35|) to (|38[) can be introduced 
into the Trueman formula (|30fl in order to express 
it in the form of Eq. © . 

The result for an arbitrary charge distribution 
follows in complete analogy to the single channel 
case discussed in Eq. (II 7[) and following, with 
the substitutions R — » (r) e etc.. The corrective 
terms are given by Eq. 118(1 with the only change 
that now a h — a 1 ^ c and b h — b^ c . The effects 
from neutral channels are negligible. 

The total level width Ti s has two components 
of comparable magnitude corresponding to the 
decay via the charge exchange and radiative chan- 
nel, respectively: 

v ls = rf s n + rj:. (39) 

They can be physically separated using the Panof- 
sky ratio P = 1.546(9) [T^]. Of special interest 
in the present context is the hadronic charge ex- 
change channel. Here the ratio b\ /a^ of Eq. (|19fl 
is bcc/a'^, since the charge exchange width de- 
pends quadratically on the amplitude l|3"2|l and 

r ls °" = y^tt = ^<t>% (0)? o Wlo (i+Sr)] 2 ■ (40) 

Here Sr is the counterpart of <5i s .i/2 of Eq. 1191) : 

Sr = b< r > +S C + \{ql + ql) h ^f + V. (41) 
z z a c,o z 

We have corrected the charge exchange amplitude 
aj? in Eq. (|40(l not only for the effective interac- 
tion energy in the charged channel (gauge term) , 
but also for the non-atomic 7r° energy in the open 
charge exchange channel l|32l) . This is justified, 
since this correction can be 'tuned' externally, for 
example by binding the proton into a potential. 
It is thus of a different nature than the non-trivial 
corrections for the mass splittings. 

3. Numerical results 

We now apply these results to the 7r~p atom. 
We assume for the moment that the correction 
for the finite range of hadronic interaction only 
enters via the range parameter b. Isospin in- 
variance is assumed for hadronic scattering pa- 
rameters appearing in correction terms. The 



e.m. expectation values appearing in Eq. l(18jl 
follow from the folded (ir~ , p) charge distribu- 
tions corresponding to the observed form fac- 
tors 19: (r) em = 0.95(1) fin; (l/r) em = 
1.48(1) fin- 1 ; (log(mr)) em = -0.687(9) with 
V c (0) = (a/r) em = 2.13(2) MeV. We use the em- 
pirical values for the range terms b^- p = b^+ n = 
-0.031(9) m- 3 ; b„- n = b^+ p = -0.058(9) m~ 3 
or b+ = -0.0044(7) m" 3 ; b~ = -0.0013(6) m~ 3 
from 7rN scattering data |20|. 

The correction terms are given in Table ^ 
For the n~p energy shift they are obtained from 
Eq. (|18J) by a two step iteration and do not require 
the knowledge of a n - p . The width corrections, 
calculated from Eqs. I|40|l and ijHJ, require that 
one knows the sign of a^ Q . We also give the cor- 
rections for the 7r + p Coulomb scattering length 
a c JH]; which is similar to the 7r _ p case, but 
for appropriate sign changes in parameters. The 
a^+p correction terms follow from our determina- 
tion of aj? c and aj? assuming them to be isospin 
invariant. The physical 7r ± n scattering lengths 
have Coulomb corrections of less than 0.1% and 
can safely be identified with the hadronic ones at 
the present level of precision. 

There is little uncertainty in any of the correc- 
tions within our assumptions. It comes mostly 
from the experimental value of the range term 
b^-p = b + + b~ . Here the b~ part contributes 
50% to the error of the energy shift and nearly 
all to that of the width. From a purely phe- 
nomcnological standpoint its theoretical origin is 
irrelevant. However, to leading order it is sim- 
ply generated by the energy dependence of the 
Weinberg- Tomozawa amplitude on the one hand 
and by the nucleon Born term of opposite sign on 
the other one (Eqs. (44-46) in Ref. 21 ), consis- 
tent with the experimental value. In the case of 
Sr, the non-atomic correction in Eq. I40JI for the 
neutral pion energy is responsible for 60% of the 
'gauge term'. 

The low-energy expansion for the K-matrix de- 
pends symmetrically on the initial and final mo- 
menta as (qf + q])/2- For the terms propor- 
tional to b~ this is explicitly the case in the non- 
relativistic limit when the initial and final pion 
are separately on the mass shell j^J • The sit- 
uation is similar for the isoscalar effective range 
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Table 1 



Coulomb corrections in percent as described in the text. S vp is included in Si s and 





Extended charge 


Renormalization 


Gauge term 


Vacuum polarization 


Total 


Sis 


-0.853(8) 


0.701(4) 


-0.95(29) 


0.48 


-0.62(29) 


Sr 


-0.427(4) 


0.701(4) 


0.50(23) 


0.24 


1.02(23) 




0.853(8) 


0.72(5) 


-1.71(29) 


0.48 


0.35(29) 



term b + . The dominant contribution is in this 
case proportional to the scalar form factor a{t) 
(see, e.g, Eq. (10.1) and following in Ref. [22]). 
The corresponding nucleon Born terms have the 
same structure. 

We have therefore a good quantitative picture 
of the precision to which the hadronic scatter- 
ing length can be extracted from the strong en- 
ergy shift in the limit of a short ranged energy 
dependent hadronic interaction. To fully exploit 
the present experimental information content, the 
theoretical corrections must controlled at least to 
0.5%. This has been achieved even using the er- 
rors of the phcnomcnological parameters. The 
true theoretical precision in our approach is far 
higher. The model with a spherical shell charge 
distribution suggests that the results are robust 
with small modifications as long as the interac- 
tion range is inside a characteristic charge scale. 

The correction terms of Table 1 applied to the 
experimental value of the pionic atom energy shift 
and width of Eq. Q give the following hadronic 
scattering lengths 

4- p = < c = 0.0870(5) m- 1 ; (42) 

<- P ^n = <o = -0.125(4) m" 1 . (43) 

Here the masses are the physical ones and the 
7r°n and decay channels are open. The re- 
sult H42fl is 1.5% smaller and outside the quoted 
uncertainty of the value 0.0883(8) deduced in 
Ref. pQ which is based on the Sigg analysis 0- We 
have made no correction in Eq. 1|42[) for the e. m. 
terms outside our Coulomb potential approach, in 
particular for the effect of the and the 7A in- 
termediate states in the direct and crossed chan- 
nels. This will be discussed in the next section. 
Our value (g3Jl for _^o n differs by 2.4% from 
that deduced in Ref. [j. 

The scattering length ajj can be analyzed 



jointly with the tt~ D scattering length to give an 
isovector scattering length (a n - p — a 7T - n )/2. We 
follow the procedure of Ref. 0] with two minor 
additions. First, the triple scattering term in the 
multiple scattering was only partly included. The 
full triple scattering term according to Refs. (2H1 
12*1) represents a contribution of +0.0027 m^ 1 (21] 
to the theoretical ir~D scattering length in 
Ref. 0] in the limit of point interactions. Fol- 
lowing the procedure in Ref. 0] we reduce the 
overall factor < 1/r 2 > D = 0.314 (25) fm -2 in 
the triple scattering term by a form factor to 
give < f{r) 2 /r 2 > D = 0.238 (24) fm" 2 . The the- 
oretical 7r~D scattering length in Ref. 0] should 
therefore be increased by +0.0019 (2) m~ . Sec- 
ond, the Fermi motion ('boost') correction now 
includes not only the dominant contribution from 
p-wave 7rN scattering, but also a smaller one from 
the energy dependence of s-wave isoscalar ampli- 
tude at threshold as first found in a chiral ap- 
proach [23] - ^ n fact, a more accurate descrip- 
tion of this correction term is obtained using the 
momentum expansion of the forward scattering 
amplitude near threshold I2£ ■ The p-wave co- 
efficient Co in Ref. 0], Eq. (14), should then be 
replaced by b + + cq with b + = (b^- p + b^- n )/2. 
The Fermi motion correction to the tt~D scat- 
tering length changes then from 0.00G1 (7) m~ 
to 0.0047(6) m" 1 . The net contribution of these 
two changes to the determination of the isovector 
scattering length is only about 0.2% of its nu- 
merical value. This is small in comparison to the 
systematic theoretical uncertainty quoted in Ta- 
ble IV of Ref. 0] . The individual contributions of 
both changes are also within this uncertainty. 

Using a n - p from Eq. (|42[l in conjunction with 
the procedure of Ref. 4 Eq.(B8), including the 
corrections just described leads to the improved 
isovector combination value {a^- n — a^- p )/\/2 = 
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—0.125 (1) m' 1 . It can therefore be directly com- 
pared to the corresponding value —0.125 (4) m" 1 
deduced in Eq. (|43|) from the charge exchange 
width rj s n . Note that when the deuteron data 
are used to extract the isovector scattering length, 
the further e. m. contributions cancel in the limit 
of charge symmetry 0]. Such is the situation 
for the leading dispersive contributions from pro- 
cesses such as the jN^A) channels and their 
crossed terms. Thus, there is no indication of 
isospin violation in the isovector amplitude at the 
present level of precision. 

4. Comparison to previous approaches 

In previous analytical approaches using wave 
functions little attention was given to the effect 
of the electromagnetic finite size effects and to 
the issue of the correct energy of the interaction. 
In a recent discussion, the Coulomb interaction 
is cut off entirely at the range of the strong in- 
teraction |Sj . Several authors starting from True- 
man |3I8| consider the influence of a threshold 
expansion for the hadronic amplitude. However, 
they incorrectly identify the momentum q c with 
the Bohr momentum, which leads to a very small 
correction of order a 2 . The correction to the scat- 
tering length proportional to (a c ) 2 is obtained to 
leading order a log a, but these previous treat- 
ments give incorrectly such terms of order a. The 
numerical consequences of this difference is small. 

As discussed in the introduction, the numeri- 
cal approach using coupled channels [5] correctly 
includes the effect of the finite size and vacuum 
polarization in the wave function correction as 
well as the renormalization term. It is, however, 
inconsistent with the low-energy expansion. In 
addition, it makes model dependent corrections 
for isospin violation and radiative decay effects 
(see also comments in Refs. The numeri- 

cal result for a single channel does not have these 
problems. It agrees with our explicit result using 
the same input parameters. 

In Refs. |9I10| . the energy shift e\ s is related to 
the scattering length in pure QCD. The calcula- 
tions are performed perturbatively to the same 
order in a as here in the framework of an ef- 
fective quantum field theory. The correction to 



Eq. (3) with this scattering length is obtained us- 
ing an expansion in in powers of a and of the 
quark masses (mj, fh u ). The leading and next-to- 
leading order terms can be given in closed form 
in terms of ttN scattering lengths [9] . The higher 
order contributions have been evaluated in the 
framework of ChPT in Refs. |9I10| . Since our 
approach is phenomenological and uses the phys- 
ical masses, the e. m. EFT corrections for mass 
splittings are implicitly included in our scattering 
length. One cannot therefore match the expan- 
sions in detail to each other. In spite of that it is 
possible to make some observations. First, the 
next-to-leading order non-analytic term a log a 
proportional to (a h ) 2 (our renormalization term) 
results quantitatively also in the EFT approach 
as in all previous potential approaches. This is 
natural, since it describes long-range physics. In- 
side our model, the charge distribution acts as a 
regulator and leads quantitatively to additional 
terms of order a. 

As discussed at the end of the previous section, 
our Coulomb scattering length A c c c in Eq. (29) 
contains additional non-potential e. m. contri- 
butions, which are an integral part of our phe- 
nomenological scattering length. Such contribu- 
tions have been calculated in a chiral quark model 
with quark wave functions |2fi| : the authors ob- 
tain expressions for the EFT empirical constants 
/i,2,3 in an approximation corresponding to jN 
and 7A appearing in intermediate states. The 
lack of knowledge of /1 is the largest uncertainty 
quoted in the EFT approach (±2.9%) pTHj . The 
authors of Ref. find that the constant f± in 
their picture is to 95% associated with the ax- 
ial baryon form factor. The corresponding direct 
and crossed contribution with as the interme- 
diate state has recently been calculated using soft 
pion techniques in the heavy baryon limit, which 
gives a 3.4(7)% attractive contribution Such 
attraction is also found in the chiral quark model 
of Ref. • This follows generically from the pre- 
dominance of intermediate states of energy higher 
than the threshold one. In Ref. the leading 
term has the symmetry structure corresponding 
to that associated with the ChPT parameter /1 . 
The next order term proportional to log 
has the same coefficient as in the corresponding 
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term 8T™ 1 of Ref . JU| , which indicates an equiv- 
alent physical expansion. 

The main difference of our approach to the 
e. m. corrections of the EFT approach |9I10) 
is the wave function and energy shift corrections 
which are linear, respectively inverse, in the e. m. 
charge radius. In spite of their ready intuitive in- 
terpretation, we have been unable to identify ex- 
plicit counterparts in the EFT approach in which 
only even powers ^r 2 ™) appear, contrary to our 
result with a leading linear power. This difference 
occurs already at the non-relativistic level [2E] ■ It 
is due to a different handling of the e. m. charge 
form factor, which is seen clearly in a configura- 
tion space representation. We modify the short 
range part of the Coulomb potential by the form 
factor and iterate its effect on the strong inter- 
action. To present order in the expansion the 
authors of Refs. |9I10| consider the effect of the 
form factor to be equivalent to an additional short 
range effective Lagrangian term proportional to 
( r em); the electromagnetic Lagrangian is other- 
wise taken to be independent of the form factor. 
This generates an e. m. contribution a/ s to the 
scattering length independent of the strong scat- 
tering, which is identical to our result to this or- 
der. However, the term in question is generated 
only by the fact that the charge distribution has 
a range. Contrary to the corresponding situation 
for a hadronic scattering length one cannot keep 
a/s (see Eq. (9)) unchanged in a strict zero range 
limit, since it would vanish. In addition, the Bohr 
wave function varies linearly with r near the ori- 
gin, contrary to the quadratic dependence for reg- 
ular potentials. Both effects make it delicate to 
iterate a finite range e. m. interaction using an 
equivalent zero range Lagrangian. We conjecture 
that this difference is only technical. The EFT 
approach should generate our terms in higher or- 
ders, since form factors are obtained using such 
descriptions. 

The 7r~p scattering length is the dominant con- 
tribution in the direct determination of the 7rNN 
coupling constant c/ 2 /(47r) via the GMO sum-rule 
as given in Eq. (4) of Ref. @|. Our result (@2J for 
the scattering length, together with the reanal- 
ysis of the -k~D scattering length as described 
above, allows an improved evaluation with a new 



value for g 2 /(47r) = 14.04(17), as compared to 
14.11 (19) found in Ref. @|. As discussed above, 
non-potential terms from the 7N (7A) channel 
contribute both to a _ and a _ ,, but cancel to 

IT p 7T a' 

leading order in in the combination required 
for the determination 4 . Additional contribu- 
tions are within the systematic uncertainties. 

5. Conclusions 

The aim of the present paper was to extract 
the phenomenological threshold scattering ampli- 
tude to high precision and model independently 
from the corresponding strong interaction energy 
shift and width in a hadronic atom such as the 
7r~p atom. Our scattering length is differently 
defined with respect to the QCD one used in 
Rcfs. |9ll()j as discussed above. We reach our goal 
within a non-relativistic picture by the following 
key observation. The Coulomb potential of the 
extended charge distribution is perfectly regular 
at short distances and it is useful to consider it 
to be an externally applied binding potential in 
addition to the hadronic interaction. This allows 
us to solve the problem exactly for a model with 
respect to which perturbations can be applied. In 
this situation the hadron masses arc the physical 
ones and the extracted hadron amplitude does 
not assume isospin invariance. The regularizing 
effect of the extended charge is an essential fea- 
ture for the understanding of the correction terms 
to the Deser-Trueman relation (J3}. It has pre- 
viously been included only in a numerical study 
using potentials 

We show then that an accurate relation can 
only be achieved if three physical effects are prop- 
erly included. First, the relevant wave function at 
the origin is not the Bohr one, but should corre- 
spond to an extended charge distribution includ- 
ing vacuum polarization. The extended charge 
distribution is at present important beyond the 
purely e. m. energy shift it produces in the atom. 
Second, the correct long-range behavior of the 
wave function induces a characteristic change of 
the wave function near the origin. The result 
of this feed-back is a quadratic correction to the 
scattering length proportional to am (a' 1 ) 2 log a 
to leading order. The leading term is known from 
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many investigations. Here we obtain a more ac- 
curate result including the terms of order a with 
little sensitivity to the model used. Thirdly, the 
low-energy expansion of the scattering amplitude 
leads to a characteristic 'gauge' correction which 
expresses that the scattering occurs at an energy 
typical of that of the extended charge Coulomb 
potential close to the origin |4I12| . This effect 
is as important as the other corrections. There- 
fore, approaches which do not respect the empiri- 
cal low-energy expansion cannot extract accurate 
values for scattering lengths from atomic data. 

The present investigation assumes that the 
charge distribution produces a Coulomb poten- 
tial that varies little over the range of the strong 
interaction. The results of our model in Eq. (|10[1 
and following suggest that the results are nearly 
unchanged as long as the hadronic range is within 
that of the scale set by the charge distribution. 
Although we have not yet investigated these as- 
pects, the results appear to be robust with com- 
pensations between contributions |25|. The sign 
and magnitude of these corrections are unlikely 
to change. Our considerations also apply to other 
hadronic systems, in particular to the 7r + 7r~ atom 
as will be discussed in a more detailed version |25) . 
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